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We determine semiclassical quasienergy spectra from peri- 
odic orbits for a system with a mixed phase space, the kicked 
top. Throughout the transition from integrability to well de- 
veloped chaos the standard error incurred for the quasiener- 
gies is a small percentage of their mean spacing. Such fine ac- 
curacy does not even require the angular momentum quantum 
number j oc 1/h to be large; it already prevails for j — 1, 2, 3. 
To open the way towards a reliable spectrum the conventional 
trace formula a la Gutzwiller has to be extended by includ- 
ing (i) ghosts, i.e., complex predecessors of bifurcating orbits 
and the possibility of Stokes transitions undergone by such 
complex orbits and (ii) collective contributions of clusters of 
periodic orbits near bifurcations. Even bifurcations of codi- 
mension higher than one must be reckoned with by accounting 
for the clusters involved through the appropriate diffraction 
integrals. 
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Since Gutzwiller's seminal work [[L||2[] it is known that 
the level density of autonomous hyperbolic systems can 
be semiclassically approximated by a sum of individual 
contributions from periodic orbits. Similarly, the spec- 
trum of integrable systems may be calculated semiclas- 
sically by EBK quantization. Gutzwiller's result was 
later extended to periodically driven systems ||j§|] whose 
stroboscopic period-to-period evolution is generated by a 
unitary Floquet operator F with unimodular eigenvalues 
e 

tr F 



lipi . The quasienergies •~p i are encoded in the traces 
n = 1,2,..., which are approximated as 
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for systems with a single classical degree of freedom. 
Each period-n orbit provides a summand determined by 
its primitive period uq, the action S, the Maslov index 
v, and the trace of the linearized map M . 

The Gutzwiller type result (|l|) can be derived from the 
integral representation 



trF" = 
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which involves, besides the Morse index p,, the action 
S(q',p;n); the latter generates the n-step map (q,p) — > 
(q',p') through S p — q, S q > = p' , where the indices on S 
denote partial derivatives. The stationary-phase approx- 
imation leading from the integral (g) to the periodic-orbit 
sum (Q) is a sensible one provided all stationary points 
are well separated. 



Generic systems, however, are neither integrable nor 
chaotic but come with a mixed phase space with stability 
islands residing in a sea of chaotic motion. The present 
Letter is devoted to the semiclassical evaluation of the 
quasienergy spectrum in that situation. In particular, 
we shall be concerned with the transition from regular 
to predominantly chaotic behavior as a suitable control 
parameter is varied. Such a transition involves complex 
changes of the stability islands when periodic orbits arise, 
disappear, or coalesce at bifurcations. The bifurcations 
generically encountered upon varying a single parameter 
and therefore said to have codimension one have been 
classified by Meyer || . The simplest type is the so-called 
tangent bifurcation at which a pair of new classical pe- 
riodic orbits is born (or coalesces and disappears for the 
opposite sense of change of the control parameter) . The 
general cases are period-m bifurcations where a "cen- 
tral" orbit of period I coalesces with satellites of m-fold 
period n = ml. At the bifurcation the n-th iterate of the 
linearized map M is the identity close to the coalescing 
orbits and gives, for one degree of freedom, the condition 
trM™ = 2. Clearly, a tangent bifurcation may be seen 
as the special case m = 1. 

As a dynamical system is driven through a sequence of 
bifurcations towards globally chaotic behavior the sim- 
ple Gutzwiller type trace formula ([j]) ceases to be a rea- 
sonable approximation to the integral (^), mostly since 
different periodic orbits approach one another so closely 
as to become incapable to yield independent additive 
stationary-phase contributions to the integral (g). Right 
at a bifurcation individual contributions to tr F n even 
diverge. To construct a "collective" contribution [^] 
in the neighborhood of a bifurcation one must approxi- 
mate the action function S in (|^) by a suitable normal 
form whose stationary points yield the cluster of classi- 
cal periodic points related to the bifurcation; the ensuing 
"diffraction catastrophe integral" then constitutes a clus- 
ter contribution to the trace tr F n in question. For some 
recent progress with several diffraction integrals relevant 
for our present study we refer the reader to |g|-|l0(| . 

When periodic orbits disappear as a control parame- 
ter passes through a critical value the nonlinear classical 
map loses a number of real solutions in favor of complex 
ones. A complex "ghost" orbit has no classical signifi- 
cance but does yield a saddle-point contribution jflj to 
the integral . In the immediate neighborhood of the 
said bifurcation the ghosts again form (part of) a cluster 
which must be treated by an appropriate diffraction in- 
tegral; then one obtains a uniform interpolation between 
the asymptotic behaviors on both sides of the bifurca- 
tion, due to the saddles for the ghosts on one side and 
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the stationary phases for the real successors on the other 
side. The simplest and best known such case arises for 
a tangent bifurcation for which the diffraction integral 
takes the familiar Airy function form ||j7]]ll|]. 

A ghost orbit often makes itself felt surprisingly far 
away from the bifurcation from which it originates, 
since the imaginary part of its action [which in princi- 
ple entails exponential suppression through the factor 
exp (— Im.S/%)] may decay slowly as one steers the dy- 
namics away from the bifurcation. We in fact find that a 
ghost frequently loses its weight through another mech- 
anism, i.e., the so-called Stokes transition after which 
the corresponding saddle of the integrand in (Q) can no 
longer be reached by deforming the original contour of 
integration to one of steepest descent without crossing 
a singularity. The transition is encountered when the 
real parts of the action are identical for a ghost (— ) 
and another "dominant" orbit in its vicinity (+), with 
Im S+ < Im S- . The phenomenon has been investigated 
in p2|p^ ] , where a uniform approximation of the suppres- 
sion factor is given. Incidentally, the Stokes phenomenon 
also implies that only ghosts with ImS > are relevant. 

Another surprise will be incurred below, in our search 
for reliable semiclassical approximations to the traces 
tiF n : Classically nongeneric bifurcations become quan- 
titatively important. These have codimension two, 
i.e., could be located only by controlling two parame- 
ters. Even though we shall be concerned with varying 
but a single parameter and never actually hit such a 
codimension-two bifurcation it seems typical to get suffi- 
ciently close for collective treatments of all participating 
orbits to become necessary. 

Leaving generalities for the moment we now turn to 
pursuing our goal for a periodically kicked angular mo- 
mentum often referred to as kicked top p^-p^t . The an- 
gular momentum J involved has components obeying the 
usual commutation rules [J x , J y ] = iJ z , etc. The squared 
angular momentum J 2 = j(J + 1) is conserved and with 
the quantum number j = |, 1, |,2, . . . fixed the Hilbert 
space is assigned the dimension 2j + 1. That dimension 
also plays the role of the inverse of Planck's constant such 
that classical behavior is attained in the limit j — *■ oo. We 
shall work here with the particular top whose dynamics 
is a sequence of rotations by angles pi alternating with 
torsions of strength ki as described by the Floquet oper- 
ator 
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The corresponding classical map may be obtained by 
writing out the stroboscopic dynamics of the rescaled an- 
gular momentum vector X = J/(J + 1/2) in the Heisen- 
berg picture and then, with the limit j — > oo in mind, 
degrading X to a c-number vector. The classical phase 
space is thus revealed as the unit sphere X 2 = 1. The 



classical stroboscopic map is easily written as the se- 
quence of three rotations, one about the x-axis by the 
angle p x + k x x, the second by p v about the y-axis, and 
the last one by p z + k z z about the z-axis. 
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FIG. 1. (a)-(d): Phase space portraits for the kicked top. 
The spherical phase space is parameterized by the azimuthal 
angle tp = q and the Cartesian coordinate z = p. Varying 
the control parameter k from to 5, the system undergoes 
the transition from integrability through mixed phase space 
to well developed chaos, (e): Bifurcation tree including pe- 
riodic orbits of period one (thin lines) and two (thick lines). 
Solid lines are real orbits, broken lines ghost solutions with 
ImS > 0. Ghost lines end at Stokes transitions, where ver- 
tical lines connect them to the dominant orbit. This is not 
indicated for strongly suppressed ghosts with large ImS. 

We perform a one-parameter study of the top, holding 
the pi fixed (p x = 0.3, p y — 1.0, p z — 0.8) while vary- 
ing the control parameter k = k z = I0k x in the range 
< k<10. For k — we incur a pure linear rotation, 
and thus classical integrability. Only two primitive pe- 
riodic orbits then arise, i.e., fixed points located at the 
intersections of the rotation axis with the unit sphere. 
For k — 5 the phase-space portrait in Fig. |l| displays well 
developed chaos. 

Bifurcations are found numerically by solving the equa- 
tion tr M m = 2 simultaneously with p = p', q = q' for 
the triple (q,p, k). All periodic orbits existing for a cer- 
tain value of k can be picked up by going through the 
sequence of bifurcations as k is swept up from zero to 
its current value. Figure |l](e) displays the bifurcation 
tree thus obtained, showing twenty orbits of length one 
and two subsequently used in the evaluation of tr F 2 for 
k < 10. 

The divergence of individual contributions at bifur- 
cations is illustrated in Fig. |^(a), which displays the 
quantum- mechanically exact tr F 2 {k) for j = 3 together 
with the sum of individual contributions from real pe- 
riodic orbits and ghosts, eq. (|]). Stokes transitions are 
taken into account for those ghosts that are not suffi- 
ciently suppressed by having a large imaginary part of 
the action. 

We proceed further to include collective contributions 
that regularize the behavior close to bifurcations. The 
broken line in Fig. @(b) is our result when one groups 
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the orbits according to the codimension-one bifurcations 
in which they participate (tangent bifurcations of orbits 
with primitive period one and two and period-doubling 
bifurcations of orbits with period one) and summing up 
the corresponding contributions, using the closed formu- 
lae from j|. 

A typical codimension-one cluster is that of the two 
period-one orbits that come into existence via the tangent 
bifurcation at k = 2.44. One of the orbits is unstable 
while the other, initially stable, becomes unstable in a 
period-doubling bifurcation at k = 4.30; as it does, a 
stable period-two orbit shows up as a satellite. Close 
to this bifurcation one has thus another cluster, formed 
by the satellite and the period-one orbit that changes its 
stability. Somewhere in between these two bifurcations 
one clearly has to rearrange the clusters. Unfortunately, 
the regrouping allows for ambiguities when an orbit is 
involved in several subsequent bifurcations. Whenever a 
regrouping is found necessary we choose its location along 
k so as to minimize the discontinuity in the approximated 
trace, taken as a function of k. In most cases, as in the 
example under discussion, the remaining discontinuity is 
tiny. In three situations, however, we have to avoid such 
patchwork and to enlarge our clusters to include orbits 
that are involved in subsequent bifurcations. In two cases 
the clusters stem from bifurcations of codimcnsion two, 
where one would have to control two parameters to let 
all participating orbits coalesce. 



exact 
codimension 1 
codimension 2 




FIG. 2. The real part of the quantum-mechanically exact 
tr F 2 (k) is compared with various levels of semiclassical ap- 
proximations. In (a), individual contributions from all real 
orbits and ghosts are summed up and Stokes transitions taken 
into account. In (b), collective contributions of orbit clusters 
are used to regularize the behavior close to bifurcations. Clus- 
ters connected to codimension-one bifurcations are found to 
be insufficient at k ~ 8. Enlarging the clusters further gives 
an accurate approximation. 



quently encountered for tr F 3 . It can be described by the 
normal form 



S^{q',p) = q'p-eq' - 
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with a = ±1. This describes three orbits, two of which 
bifurcate at e = 0. A uniform approximation is obtained 
by introducing into the exponent of eq. (Q) and ex- 
panding \S q ' p \ 1/2 = 1 + Aq' + Bq' 2 . Here A and B are 
determined by requiring that the resulting contribution 
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has the right stationary-phase limit as h — > 0, which gives 
three individual contributions of the type encountered 
in eq. (|l|). The integral can be expressed by Pearcey's 
function and its derivatives [ |l8| |. It turns out that this 
expression also correctly treats the Stokes transition of 
the complex saddles. 

Upon treating the codimension-two event as discussed 
we obtain the dotted line in Fig. |^(b) for the second trace, 
tr F 2 . This ultimate level of approximation reproduces 
the exact result quite nicely. 

The other codimension-two case encountered is that 
of a tripling bifurcation close to a tangent bifurcation 
of the satellite period-three orbit. This involves another 
satellite of period three that can be taken into account 
by an extended normal form, as is discussed in greater 
detail in where a uniform approximation is given. It 
becomes relevant in the evaluation of tr F 3 , as does also 
the third scenario of higher codimension where a tangent 
bifurcation of period-three orbits takes place on a broken 
torus formed by another pair of period-three orbits. The 
result for F 3 is considerably improved by treating all 
four orbits collectively, using the contribution 
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where all coefficients arc determined to yield the correct 
stationary-phase limit. 

With help of these collective contributions the traces 
tr F and tr F 3 come out with a quality comparable to 
that of trF 2 . 

In general, the Floquet operator F acts as an N x N 
In one of the codimension-two situations, endangering matrix with N = 2j + 1 whose eigenvalues e~ itpi are 



the semiclassical approximation of tr F around k : 
third orbit of equal length is found in close neighborhood 
to a pair of orbits that participate in a tangent bifurca- 
tion at k = 8.12, and the Stokes transition rendering the 
ghost orbits irrelevant occurs at k = 7.98. This type of 
clusters formed by three orbits of equal period is also fre- 



determined by the set of traces with, for integral j, 
n = 1 . . . j. For j = 3 the first three traces thus provide 
sufficient information to retrieve all seven quasienergies. 
Indeed, from these traces one obtains the first half of 
the coefficients in the secular polynomial det(,F — z) = 

(- z y 



l^n=0 a N 



using Newton's formulae |19|; the 
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second half follows from the unitarity of F which entails 
the so-called self-inversiveness [^0|,£l], ajv_ n = a*ajv. 
We benefit from the fact that ajv = detF, needed to 
exploit the self-inversiveness of the polynomial, is acces- 
sible semiclassically: For the top, det F factorizes into a 
product of determinants of pure rotations and torsions, 
and each integrable factor can be treated individually 
by EBK quantization, which even gives the exact result, 
detF = exp[-U(j + l)(k x + k z )]. 
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FIG. 3. (a) Exact and semiclassical level curves as a func- 



tion of k. (b) Relative error 



2lr/(2j+l) 



of the quasienergies. 



The quasienergies ifi(k) from the semiclassically ap- 
proximated secular polynomial are plotted together with 
the exact ones in Fig. ||(a). Noteworthy is the coalescence 
of two semiclassical phases mimicking a close avoided 
crossing of exact ones. The corresponding eigenvalues 
cease to be unimodular there. Self-inverse polynomials 
are capable of such behavior; to impose unitarity on them 
certain additional restrictions have to be fulfilled by the 
set of traces which are not automatically respected by 
the semiclassically approximated ones. As a quantitative 
measure of accuracy we employ the standard deviation 
between exact and semiclassical quasienergies, 
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which is shown in Fig. |^(b). The error is a small single- 
digit percentage of the mean spacing 2ir/ (2j + 1), a little 
higher only in the short k interval where the two semi- 
classical phases mentioned are degenerate. 

Similarly small is the error we have incurred for j = 2 
and even j = 1. Semiclassical behavior of the spectrum 
begins to prevail for surprisingly small values of j indeed, 
fortunately so before with increasing j the infamous ex- 
ponential proliferation of periodic orbits would render 
semiclassical work cumbersome. 

Our accuracy is comparable to the one previously 
found with different semiclassical strategies that avoid 
periodic orbits fLTjBg]. 



To summarize, we have demonstrated that the spec- 
trum of the kicked top, which has a mixed phase space, 
can be calculated semiclassically from periodic classical 
orbits, provided one improves in several ways upon the 
trace formula for hyperbolic systems, eq. (^). The mod- 
ified trace formula complements the contributions from 
isolated real periodic orbits by accounting for isolated 
ghosts as well as clusters of orbits associated with bifur- 
cations of codimension one, two,. . . . 
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